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Abstract 

A full one-loop calculation of neutralino S-wave annihilation into the final state 
is performed in the minimal supersymmetric extension of the Standard Model. This 
process, like the similar one with two photons in the final state, may be of importance 
for the indirect detection of supersymmetric dark matter through the very narrow 7 ray 
line that would result from neutralino annihilations in the galactic halo. 

We give the complete analytical formulas for this loop-induced process and treat 
the case of a pure Higgsino as a first application of our expressions. Predictions for the 
gamma line flux are given for the halo model which is of the form suggested by Kravtsov 
et al. and for the profile proposed by Navarro, Frenk and White. 

For heavy neutralinos, the lines from 27 and would have indistiguishable energy 
in a realistic detector, making the fluxes add and facilitating discovery. For lighter 
neutralinos, the positions and relative strengths of the two lines would give valuable 
information on the nature of the supersymmetric dark matter particles. 
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1 Introduction 



Although the hghtest supersymmetric particle, generally assumed to be a neutralino, is an 
excellent dark matter candidate it is not easy to prove the hypothesis that the massive halo 
of the Milky Way is predominantly composed of such particles. Due to the weak interaction 
strength, direct detection of neutralino scattering requires detectors of a sensitivity which 
only recently has started to come close to that needed according to theoretical predictions. 

In fact, even if there would appear a direct detection signal, it would most likely take a long 
time before it is established as a discovery due to the rather featureless recoil spectrum and 
the weak temporal modulations caused by the Earth's motion in the "wind" of neutralinos. 

A similar lack of distinctive signature is plaguing indirect detection of antiprotons, positrons 
and continuum gamma emission from neutralino annihilation in the galactic halo. A better 
signature is provided by high-energy neutrinos from neutralino annihilation in the centre of 
the Sun or the Earth, but km^-scale neutrino detectors may be needed to scan a substantial 
portion of supersymmetric parameter space.^ 

An excellent signature is given by a gamma ray line from neutralino annihilation into a 
two-body final state with at least one photon. Since the neutralinos annihilate with galactic 
velocities, i.e. non-relativistically, the photon in the process X + X ^ 7 + (where x is the 
neutralino and A"" is a neutral particle) will be emitted with energy 

E.=Ar-^. (1) 

In a recent paper we have performed a full one-loop calculation for the case ~ 7, i.e. 
Xi + Xi ^ 7 + 7: where Xi is the lightest of the four neutralinos in MSSM, the minimal 
supersymmetric extension of the Standard Model of particle physics. There it was shown that 
the 77 process can be at least an order of magnitude larger than earlier simpler estimates 
indicated, especially for the very heavy, pure higgsino case (however, even this larger rate is 
several orders of magnitude below that needed to explain the possible structure in the multi- 
TeV gamma spectrum recently reported in unless the halo dark matter is very clumpy 

i). 

As there soon will exist new gamma-ray detectors, both space-borne and on Earth, with an 
order of magnitude larger detection area than that of existing ones, it is of importance to 
investigate other neutralino annihilation channels which give a nearly monoencrgetic 7. In 
the MSSM, the only other such final state is Z^j Q. Note that since the annihilation in the 
halo takes place in the S partial wave, helicity conservation forbids the -ff°7 final state with 
i = 1, 2, 3 representing the three neutral Higgs bosons. There may be non-perturbative final 
states, such as Vj with V a vector meson made from a qq pair, but these generally have a 
very small branching ratio ||, ^ . 

As we will see in the case of a pure higgsino, the Z^j and 77 cross sections are of about equal 
magnitude. This gives an interesting way to verify the gamma ray detection of neutralino 
dark matter: there should be two lines with energy ratio 

Ell ^ 1 _ "^z (2) 
4M^- ^ ' 



*For a thorough discussion of the various detection methods, see jlj. 
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For a heavy higgsino, the two hnes would hardly be resolvable resulting in a further increase 
of the line strength compared to the results of For a lower neutralino mass (below a 
few hundred GeV) the two lines would be resolved. The relative strengths of the two lines 
could give a handle on the composition of the neutralino. This comes about because despite 
the fact that the two processes are closely related there are some differences (e.g. the non- 
diagonal couplings of Z*^ to charginos and squarks) which depend on composition. We leave 
the detailed investigation of this for future work, however. 



2 Cross section 

The process of neutralino annihilation into a photon and a Z*^ boson 

xUx"i^l + Z° (3) 

was briefly discussed in Q] but is examined here for the first time in a full one- loop calculation. 
We focus on the case of massive non-relativistic particles in the initial state, as this is the 
appropriate limit for neutralinos annihilating in the galactic halo (the predicted value of the 
velocity for dark matter particles in the halo is of the order of 10""^ c). The outgoing photons 
are nearly monochromatic, with energy 

K = _ (4) 

The steps we follow to compute the cross section are essentially the same as those described 
for the 2 7 case. As the P wave is highly suppressed, it is a very good approximation 
to assume that the neutralino pair is in a "'^Sq state and to calculate the amplitude using 
the corresponding projector [Q. We notice that it is convenient to choose the non-linear 
gauge introduced in Ref |^ which has the peculiarity of having a vanishing coupling for the 
J vertex (G^ is the unphysical Higgs or charged Goldstone boson). We reduce in 
this way the number of diagrams involving and loops. The diagrams giving a non- 
vanishing contribution to the cross section at the one loop level are shown in Figs. 1-4. Note 
that for a generic supersymmetric model the vertices fi fj and Z'^ xt Xj ^-re non diagonal 
respectively in the sfermion and the chargino mass eigenstates. 

The computation of the loop diagrams is greatly simplified after realising that the four-point 
functions which arise from box diagrams can be rewritten as linear combinations of three- 
point functions. This procedure was already exploited in [|| and can be applied regardless 
of the final state, because the particles in the initial state, being identical particles at rest, 
have equal four-momenta. It is easy to verify that in this case it is possible to find a linear 
combination of three or four factors in the denominator of each four-point function which 
is independent of the momentum flowing in the loop and formally reduce the calculation to 
integrals of three-point functions. 

We keep, where possible, the notation introduced in Ref. The amplitude of the process 
is factorized as 

6 ~ 

-4 = ^-^=— J e (£1,62,^1,^2) -4 (5) 
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where ei and ki are the polarization tensor and the momentum of the photon in the final 
state, whereas £2 and k2 refer to the polarization tensor and the momentum of the Z'^. The 
cross section multiplied by v, the relative velocity of the neutralino pair, is given in terms of 
A by the formula 



vaz-, 



32 7r4 



M4 



A 



(6) 



We have identified four classes of diagrams. The total amplitude is given as: 
A = Aff + Ah+ +Aw + Ag- 

For each contribution we have separated real and imaginary parts. The results are not as 
compact as those obtained for the process of neutralino annihilation into two photons. This is 
due mainly to the presence in each class of diagrams of four mass parameters. The complexity 
of the calculation, and the many parameters entering the supersymmetric models means that 
it may be difficult to check our results numerically in an independent calculation. Therefore, 
we want to give all the analytical expressions despite the rather lengthy formulas. We have, 
however, chosen a notation that facilitates a comparison with the previous 27 results B in 
the limit when Mz 0, and only diagonal terms are kept in the fi fj and Z^ xT 
vertices. 

1) Contribution of the fermion-sfermion loop diagrams. 

This class of diagrams is shown in Fig. 1. The sum over / includes the quarks and the charged 
leptons, the sum over fi and fj the corresponding sfermion mass eigenstates. 



ReAfj = Y^ 



- m| /4 



/ ^ 



E 



El [a, b, c/4) + E2 ii (a, b, c/4) 



- El li' (a, 6, c/4) + {El +E3 + Ea) /g^ (a, 6, c/4) - ^ (a, 6, c) 



E 



{E^ +Ee+ E7) (a, 6, d, c/4) + (^5 +Eq- E^) 1^% (a, b, d, c/4) 



ff t2, 



+ Es V (a. b, c/4) + Es h (a, 6, d, c/4) + 1/ (a, 6, d, c) 



E 



- M2 - m|/4 



EgIi^a,b,c/4) 



(7) 



ImA 



if Ml -ml /A 



[Y.Ei+EAj2{b,c) + J2(iEi+E3+E^)J3{a,b,c) 
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+ ■ 



C/4: 



(l-4m^/m|) 



(8) 



Cf is tlie cliargc of the fcrmion in units of the electron charge (— e), Cf is the color factor equal 
to 3 for quarks and to 1 for leptons. The functions / and / which appear in the real part 
arise from the loop integrations and the J functions are the corresponding contributions in 
the imaginary part; their explicit form is given in Appendix A. The coefficients E are listed 
below: 



E, = - 



E3 



4(1 + a- 6) 



Eo 



, qb 



4 1 - 6 + c/4 
1 {h-c/A)SK 

4 1-6 + c/4 



£4 = - 



c 



{l-b + c/A)S 



ff 



16 (a - c/4)2 

// 



1 (l/2 + d/2-c/4)5;^ 

Ek = — - 



Eq = — — 

Er = 



l/2 + rf/2-fo-c/4 
c (l/2 + rf/2-&-c/4)5;^- 



32 (a - c/4)2 

1 {i-d)s;^. 



Eg = 

Eg = 



16 a - c/4 

// 



1 bS^ 



8 l/2 + d/2-6-c/4 



4 m| 4 (4 M2 - TO^o ) 
where we have defined: 



« = TTT- 

fi 


6 = 


9 
f 


^_ < 


d = 


Mj 



fi fi 



^ff = \ (4/1 9f:n + 4n Sfji) ■ iszff + 9zn) 

qa ^ ^Zfifj / L L* , R R* \ 
ff 2 Vfjfi^fifi^fjfi^fifiJ 
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^ff - 2 (^/./i ^/'/i ^ ^^/i ^A/i ^^/■^■) 
^ff = 2 fe/i ^^^^ ^ ^f-n 4*n ^^ff) 

Gzf = 2 {Qzii - 9zii) {pzff - 9zff ) 
Gn^'f = 2 (^^^Jii ~ ^H°ii) i^Hlff ~ 9H0ff) {dzff + 9zff) 
The index 1 is referred to Xi- 

2) Contribution of the chargino-Higgs boson loop diagrams. 

This class of diagrams is shown in Fig. 2. The sums over xf ^^'^ Xj involves the two chargino 
mass eigcnstates. 



ReAH+ = 



^ M2 -m|/4 



E 



Fi li (a, h, c/A) + (F2 + Fs) 1^ (a, b, c/4) 
{Fi + Fw)hHa,d,l,c/A) 



+(F5 + Fn) (a, 1, rf, c/4) + F4 (a, b, c/4) 

+ F5 (a, 6, d, c/4) + (F5 + Fe + F7 + Fg) 72^1 (a, 6, rf, c/4) 

+ (F4 - Fe + F7 + Fg) (a, 6, d, c/4) + Fg Ig^ (a, 5, c/4) 
+ F9/34(a,6,rf,c/4) - ^7/(a,6,d,c) 



(9) 



The coefficients F are listed below: 



^ Al-b + c/A 

p ^ l {h-c/A)S^^ 
^ 4 1-6 + c/4 

c (l-6 + c/4)g^"„ 

^ 16 (a - c/4)2 

^ 8 l/2 + d/2-o-6 

^ ~ 8 1/2 + d/2 - o - 6 

^ 8 a -c/4 16 a -c/4 

c (l/2 + d/2-6-c/4)g^V 
^32 (a - c/4)2 



5 



_ 1 (l/2 + rf/2-c/4)^ 
* 8 (l/2 + rf/2-6-c/4) 

F = 1 ^ 

® 8 (l/2 + d/2-6-c/4) 

^ _ 1 VdJz + ^z 1 V^^go 
^° 16 c 8 4a-m?,o/Af2+ 

~ 16 c 8 4a-m^o/M2+ 

and we have defined: 



a = 



c = 



X? 

xj 


6 = 


xr 


2 


d = 




xr 


xr 



•^xff ~ ^ (s'h+h ^H+ii + 9H+ii) 

^XH — 2 \^H+lj 9h+1i 9z3i + 9H+lj 9H+li 9zj 

^xH = 2 (fff+ij 9H+ii 9zji + 9H+ij 9H+ii 9zji) 

^xH = 2 (^^i^+ii 9h*u 9zji + 9H+ij 9h+u 9zjij 
Sz = (sill - fff ii) {9zij 9zji - 9zij 9zji) 

Dz = (sill - 5fll) {9Z^J 9Z,^ - 9Z^, 9z,^) 

a _ ( R L \ ( L L R R \ 

^H" - \9h'^ii - Sffgii J \9Hlij 9z]i - 9H°ij 9zji j 
Dho = (9HO11 - 9hoii) (9H«ij 9zji - 9H«ij 9zji 



where the indices i and j are referred to xf xf 

3) Contribution of the chargino-W boson loop diagrams. 

This class of diagrams is shown in Fig. 3. 



Re Aw 



23 

xt 



to|/4 



Gi/i3(a,6,c/4) + G2/2'(a,6,c/4) 
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+ Gs /g^ (a, b, c/4) + (G3 + G4 + Gg) (a, 6, c/4) 



C o 



/4 (a,6,c) 



E 



G5 /i^ (a, 1, c/4) + Ge (a, 1, rf, c/4) + G5 4^ (a, 6, c/4) 



+ Ge V (a, 6, d, c/4) + (Gg + G7 + Gio + Gn ) 1^% {a, 5, d, c/4) 
+ (Gs + G7 + Gio - Gn) 72^2 (a, 6, d, c/4) + Gg Ig^ (a, b, c/4) 

+ Gs [a, b, d, c/4) - li {a, b, d, c) 



ImAw = -^Yl Ml - m|/4 ^(l-^^/^x) 



The coefficients G are given by: 



Gi 
G2 
G3 
G4 



Ge 
G7 
Gg 

Gg 



1 + a - 5 

1 3 5^V-4V^ jJ^V 

2 1 - 6 + c/4 
(l-g + b) 

1 + a - 6 

1 (2 + 6 - c/4) g^V - 4^5 D^V 

2 1 - 6 + c/4 



1 2v^D^V - (Vd/2 + 1/2) 



2 l/2 + rf/2-a-6 



6 

xW 



Gio — — 



1 2 i?^V - (\/^/2 + rf/2) 5. 

2 l/2 + d/2- a-b 

ly/E{^/d+l) - iVd + d/4: + 1/4 + c/8) S^V 

~ 2 l/2 + d/2- 6 - c/4 

1 (y/d + 1) g^V - (Vrf + b/2 + c/4) 

2 l/2 + d/2-6-c/4 
c (1-6 + c/4) g^V 

8 (a - c/4)^ 
c (1/2 + (1/2 - b - r/4) S 



h 



16 (a - c/4)' 

Gn = - 



1 (1 - d) 



(a -c/4) 



where we have defined: 



X? 


h = 








m| 


d = 






xj 



^xW ~ ~ 2 {Swu 9wii + 9wii 9wii) 
_ gcosew ,L , „K „L*\ 

^xW ~ 4 (fi'w'ij fi'mi + 9w\j 9w\i) ■ (ffiji + flfji) 



4) Contribution of the chargino-unphysical Higgs loop diagrams. 

This class of diagrams is shown in Fig. 4. 



ReAc = X! 



1 



-M2-m|/4 



(ffi + iJio) (a, 6, c/4) + (i/a + i/a - i^io) /g' (a, c/4) 



+ ^//(a,6,c) 



Ix (a. c^, 1, c/4) + (a, 1, c/A) 



+ Hi li (a, 6, c/4) + ifg I2 (a, c/4) + {H^ + He + Hr + H^) /sfi {a, b, d, c/4) 
+ {H4 -H6 + Hr + Hs) k% (a, 6, rf, c/4) + ifg li [a, b, c/4) 



5'& 

+ HsI,^{a,b,d,c/A) - -2f I^^a,b,d,c) 



(12) 



For i £ {1--9} the coefficient Hi is obtained from the corresponding substituting S^^j and 
D^jj respectively with S^q and i^^o- Hio is given by: 



4 1 - 6 + c/4 



a,b,c,d are defined as in the case of the chargino-W boson contribution and we have intro- 
duced: 



"^XG — — 2 — \3Gli 9gu + 9gu 9Gli) 
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1 
2 




yGij .9Gii yzji j 


^XG - 


2 




' .9Gii yGii yzji j 


q C _ 


1 

2 




9Glj 9Gli 9Zji) 


DxG = 


1 
2 


(ffClj 9Gli 9zji ^ 


' 9Glj 9Gli 9zji) 


^XG - 


9 
2 


(swii ffcii + 9wii 9gii) 


^XG - 


9 
2 


(SgIi ffVFli + 9Gli 



The coupling constants for left and right chiral states, and g^, are written in the conven- 
tions adopted in the PhD Thesis of Edsjo ||]; all of them are defined therein but gQ^^ and 
g^ii which are given in and 

gZGG = 9ZH+H+ 

Note that a sign factor is implied wherever a square root of a mass parameter appears. For 
instance, in the third class of diagrams 

= sign (m^o /M^+ ) • (13) 

In this way we take into account the fact that, if the convention of real neutralino and 
chargino mixing matrices is adopted, the mass eigenvalues can be either positive or negative. 

We have checked the result of the imaginary part applying the Cutkosky rules. These provide 
a general scheme to extract the imaginary part which is the contribution to the amplitude 
deriving from physical intermediate states in the loop diagrams (for a review see for instance 
Ref ) . It is sufficient to examine if in each diagram it is possible to identify two separate 
processes that are both kinematically allowed. 



3 Pure Higgsino limit 

The calculation of the annihilation rate involves many subtleties. Although all integral ex- 
pressions given here are convergent, they often individually have large values with severe 
cancellations between different terms. We have checked that our results reproduce the cor- 
responding 27 ones in the limit Mz — > 0. 

As there has recently been some interest in the TeV region due to a claim in the literature for 
a possible structure in the gamma ray energy spectrum we give as a first application of 
our results a treatment of the pure higgsino case[], leaving a more complete phenomenological 
analysis for a future work. 

We now present the results for a pure higgsino initial state. Apart from small contributions 
from the fermion diagrams, only the diagrams in Fig. 3 give an important contribution. The 

t The demand of not overdosing the Universe generally means that such a heavy neutralino has to have a 
very large higgsino fraction . 
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mass of the heaviest chargino goes to infinity, so the indices i and j refer to the lightest 
chargino. The latter is nearly degenerate in mass with Xi- We find a compact expression 
which can be easily implemented numerically. 



ReA = 



1 



to|/4 



Ci (a, 6, c/4) + 2^5/1^(0,1,2/4) 



+ (C2 + 2 C5) li {a, b, c/4) + 2 (C5 + Ct + do) 4' («, b, c/4) 
+ (C3 + 2 Cs) (a, b, c/4) + (C3 + C4 + Cg) /g^ (a, 6, c/4) 



C a 

^ h (a. c) 



06 



(a, 6, 1, c) 



(14) 



/m^ = -^— 2 ^CiJi(a,5) e{\-m'^/Ml) (15) 

The coefiicients Ci are obtained from the corresponding Gi with d=l. The degeneration in 
mass between the neutralino and lightest chargino gives a ~ 1, while the couplings are: 



9' 

4 COS 9w \ 2 



In Fig. 5, we show the values of vaz-y obtained from these expressions, as a function of 
higgsino mass. As can be seen, the maximum value is around 3.6 • lO^^'* cm'^s^^ for a 
higgsino mass around 140 GeV. For higher masses, the value of vaz-y reaches a plateau of 
around 0.6 • 10^^* cm'^s^^. This interesting effect of a non-diminishing cross section with 
higgsino mass (which is due to a contribution to the real part of the amplitude coming from 
diagrams 3e and 3f) was found also for the 27 final state in the corresponding limit, which 
gave a value of 1 • 10~^^ cm'^s^^ ||^. In the same Figure, the contribution from the imaginary 
part is also shown. This drops quite fast for the highest masses, and constitutes a unitarity 
lower bound on the cross section which agrees with the one found in H . 



4 Halo models 

An uncertain element in the calculation of the absolute gamma line fiux from neutralino 
annihilation in the galactic halo is the detailed dark matter density distribution. Since the 
integration along the line of sight involves the square of the density, any density enhancement 
may affect the predicted flux appreciably. 

Available observational data give only very poor constraints on the mass distribution within 
the Milky Way jl^. Assuming that dark matter profiles are of a universal functional form, 
we will examine two profiles derived with N-body simulations of hierarchical clustering in 
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cold dark matter cosmologies and fitted to a sample of dark matter dominated dwarf and 
low-surface brightness galaxies. We consider among the general profile family 



^^'^^ ^ {r/ap [1 + (r/a)"]('3-7)/« ' ^^^^ 

the Kravtsov et al. profile defined by (a,/3,7) = (2,3,0.2) with a very mild singularity 
at the galactic centre, and the Navarro et al. profile which has {a,(3,j) = (1,3,1). 
There are models |l^ which have a more singular behaviour near the galactic centre, and 
which would give enormously enhanced rate in that direction. However, there is observational 
evidence against such singularities from cluster gravitational lensing and the rotation curves 
of dwarf spiral galaxies ||l^. On the other hand, the discrepancy between the 1/r central cusp 
of the Navarro et al. profile and the experimental data from the dwarf spheroidal DDO 154 
has been explained in Ref. [ pT| assuming an additional component of dark baryons. 

We fix the normalization of the two profiles assuming that the dark matter density at our 
galactocentric distance, R ~ 8.5 kpc, is po ~ O.SGeV/cm"^, and we choose for the Navarro et 
al. profile a ~ 25 kpc (appropriate value for the Milky Way in an = 1 cosmology ||l^), while 
for the Kravtsov et al. we fix a ~ 11 kpc. We consider the profiles valid up to the capture 
radius of the black hole at the galactic centre ( ^ 0.01 pc for a mass of Sgr A* M ~ 10® M0). 

The gamma ray flux from the Z7 process is given by 

2 



X line of sight 



1 87 10-" ( ^'^^^ ^ /lOGcVy / R \ 
V 10-29 cm3s-i/ V ^4 / V8-5kpcy 



X 

2 

^" J(V') cm-^s-isr-i (17) 



0.3GeV/cm' 



where ip is the angle between the direction of the galactic centre and that of observation and 
J {ip) is a dimensionless function defined as 



JW = ir7^[, P\l)dm. (18) 

' /ine of sight 



R Po 



We have computed the function J{tp) numerically and the result for the two halo models 
considered are given in Fig. 6. As can be seen, for angles larger than 60 degrees the two 
curves nearly coincide and this is true even if we vary the parameters around the chosen 
values. Towards the galactic centre the Navarro et al. profile gives a huge enhancement of 
J; these values should be compared to the value 5.7, which is the maximal value (in the 
= direction) of the corresponding function for a halo described by an isothermal sphere, 
{a,f3,j) = (2,2,0), with a = R ||l^. For the Kravtsov et al. profile, the maximal value is 
J(~ 0) ~ 10. 

We are now in the position of being able to give a firm prediction for the total line flux of 
7 rays for TeV-scale higgsinos. Since any presently conceivable detector will have an energy 
resolution not better than one per cent or so, one has to add the contributions from the 27 
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annihilation calculated in and the Z7 process calculated here for the first time. With the 
Kravtsov et al. density profile, the flux in the direction of the galactic centre is given by 



$^°*(0) 



cm3s-ij (v^y 
4.10-ri^Vcm-VW-. 



— ? —1 —1 
cm s sr 



(19) 



Although this flux is much higher than previous estimates have indicated, it is still far from 
what is needed to explain the structure claimed in |^ . 

The new generation of large area air Cherenkov telescopes may have a sensitivity that comes 
close to the flux predicted in Eq. Moreover, exploiting their very small angular ac- 

ceptance (~ 10^'^sr), they will be able to investigate on the possible large central galactic 
density enhancement as predicted by the Navarro et al. profile, which could dramatically 
increase the supersymmetric dark matter discovery potential of these telescopes. 



5 Conclusions 

The results for the Z7 process show many similarities with the corresponding 27 results. For 
heavy neutralinos, the two lines would not be resolved so the line strengths would add. For 
neutralino masses below a few hundred GeV, two distinct gamma lines could be detected 
permitting a "spectroscopy" which could give valuable information on the nature (mass and 
composition) of the dark matter particles. The analytical formulas given here should be of 
great help to extract this information, if these gamma lines are detected. 

Future gamma ray detectors may be sensitive enough to discover supersymmetric dark matter 
for favourable values of supersymmetric and halo parameters, although the uncertainty of 
the latter may make exclusion of dark matter models from the non-observation of a signal 
difficult. 
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Appendix A 

In this Appendix, we define the functions needed to give the expression for the cross section. 
We start by defining the auxiliary functions 

slog{ri,r2,r3;x) = log [|-ria;^ + (ri + r2 - r^) x + raj] (A.l) 
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and 



flog{ri,r2,r3,r4,r5;x) 



X - [(7-3 + r^) /2 - r2 - r^] / (ri + r^) 
■ {log [j-rix^ + (ri - r2 + ra) X + r2|] 

- log [\r^x^ + 1/2 (rg - a; + 1/2 (rg + r^) - r^\]] . (A.2) 

These enter the following integrals, which cannot be performed analytically in terms of ele- 
mentary functions. We have chosen the integral expressions rather than rewriting everything 
in terms of dilogarithms, because the expressions are more compact and possible to integrate 
directly numerically. 

/•I (^2; 

Ii {ri,r2,r:i,r4) = / [s^og (-4ri, r2, rg; a;) - sZog (-4 r4, r2, rg; x)] (A. 3) 

Jo X 

Ii{ri,r2,r3) = (ri, r2, r2, rg) (A.4) 
I,^ (n, r2) = (ri, r2, r2, 0) = h (^i, r2) (A.5) 

{ri,r2,r3,r4) ^ / — [slog {n - 2r4,r2,r3;x) - slog {-ri,r2,rs;x)] (A. 6) 



/2'('^l:'^2,'^3) =/2'(^i>'^2,l,r-3) (A.7) 
I2'' (ri, r2) - l2^ (n, r2, 1, 0) = h (ri, ra) (A.8) 

h'^ {ri,r2,r3,r4) ^ / — [s^og (ri - 2 r4, rg, x) - s^og (-ri, r-g, r2; x)] (A. 9) 

Jo X 

hHn,r2,r3)^I:,^in,r2,l,r3) (A.IO) 
/g^ (n, r2) = /g4 (n, r2, 1, 0) = /g (n, r2) . (A. 11) 

11 iri,r2), I2 {ri,r2) and /g (ri,r2) are the functions defined in Appendix A in Ref.|^. The 
notation here and in the following has been chosen such that when Mz 0, the integrals 
have a simple relation to the ones in For example, the 1^ function below also reduces to 

12 or /g in the appropriate limits. 

(ri,r2,rg,r4,r5) = / dx [ flog {n - 2r5,r2,r3,n,r5; x) 
Jo 

~flog{-ri,r2,r4,r3,r5;x)] (A. 12) 

ha iri,r2,r3,r4) ^ I {ri,r2,r3,l,r4) (A.13) 



13 



^2^2 {ri,r2,r3,r4) = 7 ^ (n, r2, 1, n, n) 
/2^(ri,r2,r3) =/^(ri,r2,l,l,r3) 
4^ (r-i,r2,r3) = 7^ (n, I,r2,r2,r3) 



(A.14) 
(A.15) 
(A.16) 



For the pieces that can be integrated analytically, it is convenient to introduce the function 



In 



K (r, A) 



l+(r/7A) 
l-(r/VA) 



if A > 



-A arctan (:^7rs) if A < 



(A.17) 



(a, 6, c) 



1 



a_c/4 iK{l + a-b-c/2,Ai)- K{1 

-K{1 +a~b,A2)+K{l-a- b, A2) 
+K{c, c(c - 4 6)) + (1 - 6 + c/4) log 6] 



6 + c/2,Ai) 



(A.18) 



I4 (a, 6, d, c) = 



K(l + a - 6 - c/2, Ai) + K{d + a-b- c/2, Ai) 



2 (a - c/4) 

-if (1 - a - 6 + c/2, Ai) - K{d -a-b + c/2, Ai) 
-K{1 + a- b, A2) - K{d + a ~ b, A2) + K{l-a-b, A2) 
- a - 6, A2) + i^Cc + 1 - d, A3) +K{c-l + d, A3) 
+(1/2 + (i/2 - 6 - c/4) (2 log 6 - log d ) ] 



where we have defined 



Ai 


= (a - c/2 + 6 - 


lf+Aia- 


c/2) 


Ai 


= (a - c/2 + 6 - 


df + 4 (a - 


c/2)d 


A2 


= {b-a-lf- 


4a 




A2 


= (b~a-df- 


Aad 




A3 


= {l+d-cf- 


4d 





Finally, for the imaginary parts we need the following functions: 



Ji (a, b) = log 



J2 (6, c) = log 



1 + Vl-5/a 



1 + v/l-46/c 
1- v'l-46/c. 



(A.19) 



(A.20) 



(A.21) 



(A.22) 
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Js (a, b, c) = log 



(l-b + c/4:+{a- c/4) ^1 -4Vc 



1 - 6 + c/4 - (a - c/4) ^1 -4 6/c 



(A.23) 
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Figure 1: Fcynman diagrams included in the computation oi Aj^. Diagrams with exchanged 
initial states are not shown. 




Figure 2: Feynman diagrams included in the computation of >liif+ . Diagrams with exchanged 
initial states are not shown. 
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Figure 3: Fcynman diagrams included in the computation of • Diagrams with exchanged 
initial states are not shown. 




Figure 4: Feynman diagrams included in the computation of Ac- Diagrams with exchanged 
initial states are not shown. 
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-27 




Higgsino mass (GeV) 



Figure 5: Annihilation rate of pure higgsinos into a photon and a T? boson obtained in this 
work (solid line). Also shown is the luiitary bound coming from keeping the imaginary part 
only, which agrees with the result of Ref. [4] (dashed line). 
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Figure 6: Function J{ip) as defined in Eq. (18). In the detail the discontinuity in the two 
curves is due to the portion of halo which is hidden by the black hole at the galactic centre. 
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